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ABSTRACT 

The efficiency of nonmodal seif-lieating by acoustic wave perturbations is examined. 
Considering different kinds of kinematically complex velocity patterns we sliow that 
nonmodal instabilities arising in these inhomogeneous flows may lead to significant am- 
plification of acoustic waves. Subsequently, the presence of viscous dissipation damps 
these amplified waves and causes the energy transfer back to the background flow in 
the form of heat; viz. closes the "self-heating" cycle and contributes to the net heating 
of the flow patterns and the chromospheric network as a whole. The acoustic self- 
heating depends only on the presence of kinematically complex flows and dissipation. 
It is argued that together with other mechanisms of nonlinear nature the self-heating 
may be a probable additinal mechanism of nonmagnetic chromospheric heating in the 
Sun and other solar-type stars with slow rotation and extended convective regions. 
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1 INTRODUCTION 

It is widely known that stars with outer convective zones, in 
terms of their chromospheric line emission, exhibit twofold 
behavior. Younger, rapidly rotating stars have the power of 
the chromospheric emission higher by 2-3 orders o f magni- 
tude than slowly rotating, older stars (|Noves et al.|[T984i ). It 
is believed that the emission in these two classes of stars is 
related with different physical processes: the low emission 
of slow rotators represents the effect of acoustic, nonmag- 
netic heating, while the higher emiss ion of fast rotators i s 
caused by effects of magnetic activity (jSchriiver et ahlFigSgl ). 
In the last decade it became increasingly clear that var- 
ious slowly rotating, solar-type stars, with well-developed 
convection zones, have similar chromosphe ric emission. The 
available observational evidence suggests (|Goodmanl |2004| ) 
that this emission is caused, at least partially, by the upward 
propagation and dissipation of sound (acoustic) waves. 

On the other hand, statistical analysis of existing ul- 
traviolet (UV) emission data of solar-type stars across the 
H-R diagram, has shown the presence of a characteris- 
tic, well-pronounced "bas al" level of chromospheric heating 
l| Judge fc Carpenteill 19981 ) . It is commonly assumed that the 
heating may be caused by the nonlinear acoustic wave shock- 
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heating mechanism, but the dynamics of this phenomenon 
is not completely understood. Moreover, the direct exam- 
ination of Goddard High-Resolution Spectrograph data of 
several evolved stars, having similar "basal" levels of chro- 
mospheric activity, failed to detect a definite evidence in 
favor of the presence of acoustic waves, presumably a key 
component of the heating mechanism. Instead, in a number 
of counter-intuitive behavior, reminiscent of the so- 

lar transition region was detected, suggesting rather a mag- 
netic heatmg mechanism for these stars. This controversy led 
to the tentative conclusion that upward-propagating shock 
waves do not necessarily dominate the observed radiative 
losses from chromospheres of stars exhibiting typical "basal" 
behaviour. More generally, the nonmagnetic nature of the 
basal components of convective, solar-type stars has been 
called in question. 

In solar physics the problem of the chromospheric heat- 
ing has been discussed ever since it was found out that the 
chromosphere is hotter than the photosphere. Or iginally it 
was proposed (|Biermanlll946l : ISchwarzschild|[T948l ) that the 
solar chromosphere is heated via the dissipation of acoustic 
waves generated by the overshooting convection zone turbu- 
lence. A theoretical model, determining the spectrum and 
power of these waves, had been developed on the basis of the 
LighthilbStein theory l|Lighthil]||l952l : ISteinI 1967h of the con- 
vective zone acoustic wave emission. It was suggested that 
the dissipation of turbulently created acoustic waves may 
explain the observed chromospheric heating in the Sun. In 
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this model it was supposed (|Kalkofenll20"07l ) that the heating 
takes place in the granules. Al ternatively, it ma y also occur 
in the lanes between granules (jLites et al.lll99>3 ). where the 
acoustic monopole emission may induce the chromospheric 
oscillations. At any rate both these scenarios were based on 
the assumption, confirmed by various observations of the 
solar chromosphere, that the acoustic mechanism of heating 
may be important in regions where the magnetic field is of 
a little significance both dynamically and energetically. 

However, similarly to the situation with other solar- 
type stars, this explanation of the chromospheric heating 
has difficulties also for the Sun. The analysis of solar data 
from the SUMER(Solar Ultraviolet Measurements of Emit- 
ted Radiation) instrument on SOHO, analogously to the 
above-mentioned data from other slowly rotating, solar- 
type stars, revealed serious problems with the nonlinear 
acoustic-shock-wave interpretation and has underlined the 
necessity for the further investigation of this issue. Obser- 
vations with TRACE (Transition Region And Coronal Ex- 
plorer) NASA space telescope - revealed a very substantial 
(about 90%!) deficit in the energy flux required to heat the 
nonma g netic, quiet solar chromosphere ( Fossum fc Carlsso n 
l2005al lbl. I2OO6I '). Later some light on this puzzle was shed 
by numerical simulations wi th a three-dimensional hydro- 
dynamic code. It w as argued l|Wedemever-Bohm et al ] |2007l : 
ICuntz et al.ll2007l ) that this puzzling discrepancy may arise 
due to the limited spatial resolution of the TRACE, fail- 
ing to detect most of the acoustic waves traveling in quiet- 
Sun inter^^networkj^gions of the chromosphere. In particu- 
lar, [Cuntz^^eF^l] (|2007 ) discussed the problems of the limited 
sensitivity of TRACE while studying the three-dimensional 
solar chromospheric topology; they reconsidered the acous- 
tic chromospheric wave energy fiux and related this prob- 
lem with the heating and emission of chromospheric basal 
flux for T Ceti type stars. It was concluded that, contrary to 
claims by Fossum & Carlsson, high-frequency acoustic waves 
are sufficient to heat the nonmagnetic solar chromosphere. 

Recently various alternative scenarios of the chromo- 
spheric heating for solar-type, cool stars with partially ion- 
ized chromospheres have been propo sed. It was argued 
that the Farley-Buneman (|Farlevlll963l ') plasma instabilitjQ 
may provide the mechanism for conversion of the energy 
of convective m otions of neut r al atoms into the c hromo- 
sphe ric h eating f'Fontcnla"2005': ' Fontenla et al.ll2008l ). How- 
ever [Cogoboridzo ct al. (2008), taking into account the fi- 
nite magnetization of the ions and Coulomb collisions, re- 
evaluated the instability for the solar chromosphere and ar- 
gued that it could hardly be fully responsible for the ob- 
served level of quasi-steady heating. 

Therefore, the problem of the chromospheric heating 
remains an actual and unsolved problem in the physics 
of slowly rotating, solar-type stars and the Sun. Putting 
the problem in plain words, it is observationally obvious 
that the chromospheres of stars with extended convective 
zones are hotter than their photospheres. This is true for 
both younger, rapidly rotating stars with magnetically ac- 
tive chromospheres and for older stars, with slower rota- 



tion rates and hydrodynamic, nonmagnetic chromospheres. 
iKalkofenI l|2007l ) has recently considered both possible mech- 
anisms - the magnetic and the nonmagnetic (acoustic) ones - 
that might be re sponsible for the solar chromospheric heat- 
ing. He argued l|Kalkofenl [20o3 ) that the emission of the 
nonmagnetic chromosphere exhibits some characteristic fea- 
tures of acoustic waves and no signatures of magnetic waves. 

Still, in general, the problem is complicated. From the 
existing observational, theoretical and numerical studies of 
the "chromospheric heating problem" for cool, slowly rotat- 
ing, solar-type stars, it turns out that the problem doesn't 
have a single, broadly accepted solution. Rather, it seems 
likely that there may exist several different channels (mech- 
anisms) of the chromospheric heating. It is a challenging task 
to specify these mechanisms and to provide at least a quali- 
tative description of their relative importance. It seems evi- 
dent that at any rate the chromospheric heating by acoustic 
waves takes place. The question is whether nonlinear heat- 
ing by means of acoustic shock-waves is the only, or domi- 
nant, hydrodynamic channel of the heating or there are other 
acoustic processes of the energy transfer from waves to the 
medium that could lead to the chromospheric heating. 

The role of the chromospheric flows and their kinematic 
complexity is possibly one of those important issues that 
has to be tested in the chromospheric heating context. It 
is well-known that the structure of the solar chromosphere 
is approximately spherically symmetric for about the first 
2000 km. Higher layers show complex and variable fine struc- 
ture, consistin g of a wide variety of flow patterns - sp i cules, 
fibrils, surges (jBrav fc Loughhead|[l973 : ]Athavlll986l : IZirinI 
'1988!). Recent observations by Hinode strengthened this pic- 
ture , showing tha t the chromosphere ofte n hosts giant Ha 
lets iKosugi et al.l (|20o3); iNishizuka et all l|2008l ). It is now 
widely believed that the solar chromosphere, both in its ac- 
tive and quiet phasefl is highly dynamic and highly inhomo- 
geneous consisting of numerous flows ("velocity patterns") 
with different spatial/temporal time-scales and different ge- 
ometries. These flows are closely related to the chromo- 
spheric network - a web-like structure seen in the emissions 
of the red line of hydrogen (Ha) and the ultraviolet line of 
calcium (Ca II K). The network is formed and maintained 
due to the presence of bundles of magnetic field lines cease- 
lessly shuffled and reshuffled by the fluid motions. 

In particular, common ingredients of the solar chro- 
mosphere are spicules (small, jet-like ejections com- 
monly observed throughout the chr omospheric network) , 
macro- spicules and solar tornadoes (|Pike fc MasonI Il998l : 
iRogava et al]|2000l ). These flows are characterized by spa- 
tially inhomogeneous velocity fields: shear fiows with a non- 
trivial geometry and a poorly known but considerable kine- 
matic complexity. It is well-known that collective phenom- 
ena in shear flows are strongly affected by so-called non- 
modal processes, related with the non-self-adjointness of 
the o perator describing the linear dynamics of perturba- 
tions (|Trefethen et al.lll993l ). These phenomena are espe- 
cially robust and versatile in flows with nontrivial kine- 
mati cs and geo metries. In particular, it has been shown 
(Mahaian fc Rogava 1999. ) that under favourable conditions 
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a hydrodynamic, linear system may exhibit strongly pro- 
nounced unstable regimes, where due to the efficient ampli- 
fication of waves a large part of the mean flow kinetic energy 
is pumped into the unstable modes fostered by the flows. On 
the other hand, it was demonstrated that if viscous dissipa- 
tion is present, the large amplitude, nonmodally amplified 
waves undergo viscous dissipation, leading to the conversion 
of the wave en ergy into the rmal energy and corresponding 
"self-heating" (|Rogavall2004l ') of the "parent" flows. 

The physics of the self-heating in a differentially mov- 
ing medium (shear flow) is quite generic and comprises 
the following three steps: (a) The waves, originally, are 
excited spontaneously within the shear flow; (b) They 
undergo nonmodal amplification, they grow, extracting 
a part of the mean flow's kinetic ene rgy. These pro- 
cesses arc quite e fficient in neutral fluids l|Butler fc Farrelll 
[f,9,92; Chagclishvili et all [W97) and both in electrostatic 
iRoeava et al.lll997l ) and magnetized plasmas (|Rogava et al.l 
I2OO0I ). (c) Nonmodallv amplified large-amplitude waves un- 
dergo, as the final stage of their evolution, viscous decay: 
they get damped and give their energy back to the "parent" 
flow in the form of heat. 

Evidently the crucial element of this three-step process 
is the second one: nonmodal amplification of waves by the 
flow. The first and the last ones - spontaneous excitation 
of waves and their viscous damping - are routinely taking 
place in all kinds of continuous media. The possibility of the 
nonmodal amplification of collective modes in shear flows, 
however, creates a chance for the waves to amplify before 
they get damped! The energy needed for the amplification 
is drawn from the "reservoir" of the shear flow kinetic energy 
iButler fc Farrelll [1992 '). Eventual damping of these larger- 
amplitude waves and conversion of their excessive energy 
into the heat leads to the overall heating of the medium. 
Since this process converts, via the agency of nonmodally 
amplified waves, the kinetic energy of the flow into its ther- 
mal energy the process can be called self-heating. 

Self-heating was originall y described fo r acoustic waves 
in a hy drodyna mic shear flow ('RoKava"2004') , but later it was 
found (|Li et al.. .2006: Shcrgclashvih ot al. 200^) that self- 
heating occurs in magnetohydrodynamic (MHD) shear flows 
too. The particularly robust nature of t he nonmodal ampli- 
flcation in kinematically c omplex flows (|Mahaian fc Rogaval 
I1999I : iRogava et al.|[2003al ) suggested that self-heating could 
be especially well pronounced in flows with complicated ge- 
ometry and kinematics. 

Since the solar chromosphere reveals the abundant pres- 
ence of various nontrivial velocity patterns, and there is no 
reason to believe that chromospheric zones of other solar- 
type stars are any different it is reasonable to investigate 
the above mentioned mechanism in the context of the chro- 
mospheric heating for the slowly rotating, solar-type stars. 
In the present paper we consider the nonmodal self-heating 
mechanism in hydrodynamic, kinematically complex flows 
and study the energy dissipation rate of acoustic waves in 
the context of the possible relevance of this process to chro- 
mospheric flows. 

The paper is arranged in the following way: In section II, 
we derive the equations governing the linear evolution of 
perturbations of the hydrodynamic system and consider the 
self-heating occurring within a kinematically complex flow. 
In section III, we apply this model to the chromospheric 



flows and solve the set of equations numerically for different 
parameter regimes. In section IV, we discuss the obtained 
results and outline the directions of further, more detailed, 
quantitative and applicative studies. 



2 MAIN CONSIDERATION 

In order to show the efficiency of the self-heating for acoustic 
waves generated, nonmodally amplified and viscously dissi- 
pated in kinematically complex shear flows, we consider the 
following standard hydrodynamic set, consisting of the equa- 
tions of mass conservation: 

Ap + pV-V = 0, (f) 
the momentum conservation: 

DtV = --VP + i^AV, (2) 
P 

and the polytropic equation of state: 

P = Cp'+'/". (3) 

In these equations, Dt = dt + (V ■ V) denotes the con- 
vective derivative, p is the density, V is the velocity, 1/ is 
the coefficient of kinematic viscosity and n is the polytropic 
index. 

The instantaneous values of the physical variables are 
expressed as sums of their equilibrium and perturbational 



components: 

P = po + p', (4) 

V = Vo-hV', (5) 

P = Po + P'. (6) 



In the present study, for the sake of simplicity, the density 
of the fluid in the unperturbed state is assumed to be ho- 
mogeneous (po = const). 

Applying the standard procedure of linearization to 
Eqs. (|1I2I and |3]) in terms of Eqs. (|4I5I and |6]) one derives 
the following set of equations for the perturbed quantities: 

Otp' + Po(V-V') =0, (7) 

I?tV'-h(V'-V)Vo = -— Vp' + i^AV', (8) 
po 

where = a* -f (Vo ■ V) and C'i = dPo/dpo- 

Since there is no first-hand information about the veloc- 
ity fields Vo of chromospheric flows and it is impossible to 
specify 'typical' or 'standard' flow patterns we do not limit 
this study with a narrowly chosen model for a prototype 
"parent" flow. Instead we will follow the approach originally 
developed and used in lMahaian fc Rogaval (|l999l ) which ap- 
plies to a broad range of possible flow patterns. Doing so 
we expand the velocity in a Taylor series in the vicinity of a 
point A(xo,yo, zo). Preserving only linear terms we write: 
3 

V = v( A) + y (^, _ ^,0) , (9) 

1=1 

where i = 1, 2, 3 and Xj = (x, y, z ). 

In lMahaian fc Rogarol (|l999l ) it was shown that the ini- 
tial set of linearized partial differential equations (|7I8|I is 
transformed to the set of ordinary differential equations by 
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means of the usage of the following ansatz for the system's 
physical variable^ 



F{x,y,z,t)=F{t)e^'-^\ 

3 



i = l 
3 

i = l 



Ki{t)dt, 



(10) 
(11) 

(12) 



where Vi(j4)'s are the unperturbed velocity components and 
Kiitys are the wave vector components. The latter satisfy 
the following set of ordinary differential equations: 



K = 0, 



(13) 
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K 




^y,y 


Vy 


14,. 







where S'^ is the transposed shear matrix with the shear 
matrix S jMahaian fc Rogava|[l999l ) defined as: 



(14) 



with Vi^k = dVi/dxk 

Therefore, the ansatz reduces the mathematical aspect 
of the problem to the study of the initial value problem for 
the amplitudes F(t) of the spatial Fourier harmonics of the 
perturbations. 

The background flow velocity comprises stretching and 
rotation of the flow field lines in the XOY plane super- 
imposed on a spatially inhomogeneous outflow along the 
Z axis. In terms o f the shear matrix it is written as 
l|Rogava et al.ll2003al lbl): 



(15) 



Using the notation Z'^-* = df/dr for temporal deriva- 
tives, where r = tKnCs is a dimensionless time variable, 
we can write down the equations for the perturbations, i.e. 
Eqs. (|7I8|I . in the following, completely dimensionless form: 




+ k-u = 



u 



0, 

= kd- 
= 0, 



(16) 
(17) 
(18) 



where the following dimensionless notations are used: d = 
~ip'/po, u = V'/Cs, a = T./K„Cs, ai,2 = Ai.2/K„Cs, 
Ri,2 = Ci,2/K^Cs, k = K/ifn and V = uKn/Cs. Here iC„ 
is the value that will be used for the normalization of the 
wave vector. In particular, for the two-dimensional case we 
take K„ = Kx{0), while for the three-dimensional case - 
K„ = Kz{= const). 

For the sake of the forthcoming analysis it is useful to 
define the total energy of the perturbations Etot and its 
components: the kinetic, Ekin, and the compressional, Ec, 
energies: 

^ Note, that in hydrodynamics similar method, designed for the 
study of incompressible perturbations in flo ws with spatial l y uni - 
for m shearing rates, has been i ntroduced bv iLagnado et all d 19841 ) 
and lCraik fc Criminalel l ll98d) . 



u d 

Etot^Ek.^ + E^^^ + -. 

It is straightforward to verify that: 
E'-^^ = -u(S ■ u) - Vk^u. 



(19) 



(20) 



For the evaluation of the efficiency of the self-heating we 
will also use the self-heating rate, defined as (jRogava. 20041 ): 



*(r) 



k {t')u {r')dT'. 



(21) 



3 RESULTS 

From previous studies it is known that self-heating is strong 
enough only if the instabilities are robust and swift enough 
to pump energy into the acoustic waves before they are 
dumped by the viscous dissipation. At the other hand, the 
viscosity must be strong enough to ensure the feedback 
mechanism, i.e. the transfer of the energy acquired by the 
waves back to the flow via viscous damping. It is reasonable 
to expect that strongly pronounced self-heating occurs only 
when there is a cert ain balance between the ins tability and 
viscous time-scales l|Rogavall200"3 :llI et al.|[2"006h . Neither in 
the case of weak nonmodal amplification (i.e. when the waves 
do not have enough time to grow before the viscosity kicks 
in) , nor in the case of negligibly weak viscosity (i.e. when the 
dissipation is not strong enough to dump the waves and heat 
the matter), the self-heating mechanism will not be efficient 
enough to ensure a substantial heating of the flow. From 
Eqs. (|17|) . by means of a simple dimensional analysis, one 
can obtain a value of the maximum length-scale at which 
the viscous terms start becoming important: Acr ^ (iv/Ca 
(where (3 is & dimensionless parameter). Below, we will see 
that an efficient self-heating can only happen for relatively 
small A Acr length-scales. 

In this section we intend to show how strong the self- 
heating could be for flows of different geometry and kinemat- 
ics. Instead of restricting the consideration to a speciflc kind 
of chromospheric flow pattern we try to show the potential 
signiflcance of this nonmodal process for various kinds of 
kinematically complex hydrodynamic flows. 

As a flrst example, we consider the simplest, two- 
dimensional case. The Eqs. (|16m8p then reduce to: 



rf'"*"' + kxUx + kyUy = 0, 

u"^^ + aux + aiUy = kxd — vk^Ux, 
u^y^ + a2Ux — aUy — kyd — Vk^Uy, 
k^x'''' + o-kx + a2ky = 0, 
k^j^^ -\- aikx — o'ky = 0. 



(22) 
(23) 
(24) 
(25) 
(26) 



Physical parameters of chromospheric flows in solar- 
type stars may have broad ranges of variability. In the low 
solar chromosphere, for instance, the flow velocity V, lies 
in the range 10 — 25 km/s , while the temperature, T, and 
the number density, n are of the order of 5 000 — 6000 K 
and 10^^ cm~^, respectively l|Sterling|[2OO0h . The length- 
scale, R, of chromospheric structures may vary in the range 
300 — 1500 km. One can easily show that for these values 
the nonzero components of the shear matrix ~ VXx/{RCs) 
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Figure 1. Temporal behaviour of the sclf-hcating rate and the 
normahzed total energy of perturbation. The set of parameters is: 
ai = —1.3, a = a2 = 0, k^o = 1, kyo = —4, u^o = 0.2, Uyo = 0, 
T = 6000K, = leOkni, n = 10"cm-^. 



Figure 2. Temporal behaviour of the self-heating rate and the 
normalized total energy of perturbation. The set of parameters 
is: ai = a2 = —u = 1, k^o = kyo = 1, u^o = 0.2, Uyo = 0, 
T = 6000K, = 200km, n = lO" cm-3. 



[see Eq. ([29j] for Xx ^ R are of the order of unitjQ. For the 
viscosity coefficient, we use the expression that is valid for 
a weakly ionized plasma, which appears to be a good ap- 
proxi mation for the given range of temperatures (|Goodmanl 
I2OOOI ): 



(27) 



where ao is the Bohr radius and m^-the proton mass. 

It is well-known that shear flows with a consider- 
able kinematic complexity are subject to various, non- 
modal instabilities, sometimes of a parametric nature 
ijMahaian fc Roeaval ligggl ). The nonmodal evolution of 
perturbations strongly depends on the flow parameters 
and the wave chara cteristics. For instance, it was shown 
(jRoeava et al.ll2003al lbll that if F = ct^ + aia2 > 0, the sys- 
tem undergoes a linear instability, while for F ^ the flow 
remains stable and exhibits an instability only for a narrow 
range of perturbation parameters. 

For the first example, we consider the case when the 
velocity shear is one-dimensional. In Fig. ^ the time evo- 
lution of the self-heating rate and the total energy of the 
acoustic perturbation for the case of a parametric instabil- 
ity (r ^ 0) are displayed. The set of parameters in this case 
was: ai = —1.3, a — a2 = 0, kxo ~ 1, kyo — —4, Uxo ~ 0.2, 
Uyo = 0, T = 6000 K, = 160 km, n = lO" cm"^. From 
these plots it is clear that initially the unstable growth domi- 
nates over the viscous damping and the energy of the pertur- 
bation steadily increases. Gradually, however, the viscosity 
becomes more important and, as a result, the perturbation 
energy reaches its maximum level and starts decreasing until 
the perturbation completely disappears, giving back all the 
gained energy to the flow. Note that in this case the asymp- 
totic value of the self-heating rate ^'□0 = limt^oo ^(i) — 60 
is reached. 

The second example shows the case [see Fig. ©] when 
4 Here, ~ 1/Kx{0) 



the wave number vector components evolve exponentially 
(r > 0). The set of parameters is: ai = 02 = —a — 1, 
kxO = kyo = 1, Uxo = 0.2, Uyo = 0, T = 6000 K, 
Xx ~ 200 km, n = 10^^ cm~^. The qualitative behaviour of 
the system is similar to the previous case: initially the per- 
turbation rapidly grows at the expense of the background 
flow energy. But due to the exponential growth of the mod- 
ulus of the wave vector the perturbation length-scale rapidly 
decreases, leading to the inevitable domination of the viscos- 
ity, which, in its turn, eventually damps the perturbation, 
converting its excessive energy into the thermal energy of 
the flow. The figure shows that the asymptotic self-heating 
rate in this case is about 'I' 00 ~ 70. 

Observational evidence suggests that the solar chromo- 
sphere is populated by three-dimensional fluxes, viz. flow 
patterns with nontrivial morphologies. An interesting class 
of these structures a re the swirling macros picules, also called 
the solar tornadoes IPike fc MasorJ ()l998l ). characterized by 
a helical motion of the plasma. It is plausible to expect that 
similar, highly nontrivial flow patterns may exist in chro- 
mospheric layers of other slowly rotating, solar-type stars. 
These are particularly intricate examples of flows with a 
very high degree of kinematic complexity! They may host 
various kinds of str ongly-pronounced nonmodal shear insta- 
bilities [Rogw^^^ial. (2003a b). Therefore, in our analysis, 
it is worthwhile to consider three-dimensional cases too and 
to check whether self-heating processes are equally or more 
efficient in them. 

We consider the velocity configuration described by: 



V(r) = [0,rr!(r),[/(r)], 



(28) 



with f2(r) = A/r" and U = const. In this c ase the shear 
matrix can be written in the following way ijRoeava et al.l 
(|2003al )): 



S= y4(l-a)/r-" 
V Ri 

As an example, let us consider the case a 




(29) 



0.2, 



6 Andria Rogava, Zaza Osmanov and Stefaan Poedts 





Figure 3. Temporal behaviour of the self-heating rate and the 
normalized total energy of perturbation. The set of parameters 
is: a = 0.2, r = R/2, R\ = R2 = —1, fczO = —kxO = 1, kyo = —5, 
u^O = 0.1, Uyo = u^o = 0, R = SOOfcm, U = 20km/s, T = 6000/4", 
\x = 200A:m, n ■■ 



V^{R) = U, and r = R/2. Fig. ^ shows the solution of 
Eqs. (|16m8|l for the following set of parameters: Ri = R2 ~ 
— 1, kzo ~ —kxo = 1, kyo = —5, Uxo — 0.1, Uyo — u^o = 0, 
R = 500km, U = 20km/s, T = 6000if, = 200km, 
n = 10^^ cm""^. From the plots we see that the instability 
is so robust that the viscous dissipation fails to affect the 
perturbation evolution in a graphically visible way. Never- 
theless, despite the seemingly insignificant contribution of 
viscous terms, their role in the self-heating process is as 
much instrumental as in all previously considered cases. The 
figure shows that the self-heating rate at the end of the given 
time interval (r = 20) is quite high >!/ ~ 500. 



4 DISCUSSION AND CONCLUSIONS 

Slowly rotating, older stars with well-developed convection 
zones (solar-type stars) feature similar chromospheric emis- 
sion. Observations indicate that the emiss ion is related to 
the dissipation of sound (acoustic) waves (|Goodmanll2004l ') 
propagating from the photosphere to the chromosphere. 
The statistical analysis of the UV emission data from these 
stars shows the presence of a significant level of chromo- 
spheric heating (jjudge fc Carpenteilll998l ). Commonly it is 
assumed that the heating is due to the dissipation of non- 
linear acoustic waves. But the dynamics and the strength 
of the shock-heating mechanism are not well-understood 
and the nonmagnetic nature of chromospheric heating for 
solar-type stars has been recently doubted. From Goddard 
High-Resolution Spectrograph data of a number of solar- 
type stars no firm evidence was detected of the presence of 
acoustic waves. Quite on the contrary, in some cases the ob- 
servational evidence seemed to be in favour of a magnetic 
mechanism for chromospheric heating for these stars. On 
the basis of this controversy it was surmised that upward- 
propagating shock waves do not necessarily govern radiative 
losses from the chromospheres of solar-type stars. 



The situation is much the same in the context of so- 
lar physics. It is commonly believed that the solar chro- 
mosphere, in its quiet phase, in regions with negligibly 
weak magnetic field, has to be heated by acoustic waves. 
However, the energy fiux of these waves, measured in the 
upper photosphere by TRACE, has been found to be in- 
sufficient to explain the radiative e missio n fro m the chro- 
mosphere. IWedemever-Bohm et al] l|2007h and lCuntz et al.l 
(2007), employing the three-dimens ional hydrodynamical 
model l|Wedemever-Bohm et al.1l2004 ). have suggested that 
the spatial resolution of the TRACE is insufficient to resolve 
inte nsity fiuctuation s that occur on small spatial scales. 

iKalkofenI \200% has also considered this assumption, 
implying that the spatial averaging by TRACE may serve 
as a qualitative explanation for the observed acoustic fiux 
deficit. However, he found out that the standard hydrody- 
namical model was too oversimplified in the treatment of 
chromospheric energy exchange to provide a quantitative 
explanation of the suppression of t he acoustic fiuctuations. 
The m entioned hydrodynamic test (|Wedemever-Bohm et al.] 
l2007t ) was of a preliminary nature, because the hydrody- 
namic model on which the test was based involved temper- 
ature fiuctuations exceeding those of the Sun. The shape of 
the acoustic spectrum observed with TRACE seemed to sup- 
port the theory of wave generation in the solar convection 
zone even though the limited acoustic frequency range and 
the low energy flux of the observations did not allow to make 
any definitive conclusion l|Kalkofenll2007l . I2OO8I ) . That is why 
the conclusion was that "the heating mecha nism of the chro - 
mosphere thus remained an open question" | Kalkofenl|2008| ) . 

The main goal of our present study was to show that the 
nonmodal self-heating by acoustic waves may serve as an ad- 
ditional source of chromospheric heating. We demonstrated 
the possibility of the efficient heating in chromospheric fiows 
due to the presence of two basic factors: the existence of a 
spatially inhomogeneous, sufficiently complex velocity field 
(kinematically complex "flow patterns"), and the existence 
of viscous damping. 

The self-heating is efficient when some necessary condi- 
tions upon the fiow characteristics, for instance the kinetic 
energy budget of the mean flow, are met. The energy of the 
regular motion has to serve as an ample energy reservoir for 
the growing waves. Therefore, it is reasonable to consider the 
kinetic energy density of the flows in chromospheric struc- 
tures and to compare it with the energy density of unstable 
modes. 

Let us suppose that the characteristic velocity and the 
mass density of a chromospheric structure (e.g., a macro- 
spicule, a solar tornado, or a jet-like surge) are V and p, 
respectively. The kinetic energy density can be expressed 
by: 



2 



(30) 



For the initial energy of the induced mode one gets [see 
Eq. dSJ]: 



Bind = {ul + do) 



(31) 



The energy converted into the form of heat, at the other 
hand, is of the order: 
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En,at~^^{ul + dl) , (32) 

implying that each generated mode extracts a certain por- 
tion of the mean fiow kinetic energy: 

r, = ^f^{ul + dl). (33) 

For the parameters used for Fig. Q and Fig. ^ one can 
easily show that -q equals ~ 0.5 and for Fig. ~ 1, lead- 

ing to the conclusion that for the cases shown in Fig. ([1]) and 
Fig. ((2)1 about 50% of the chromospheric structures' kinetic 
energy has been extracted and given back to the flow in the 
form of heat and almost ~ 100% for the case presented in 
Fig. ([3]). It means that this process is quite efficient for the 
heating of the "parent" flow. Apparently, if a large part of 
energy is extracted from the flow patterns they could easily 
be destroyed by the self-inflicted process of self-heating. This 
circumstance, in turn, could contribute to the statistically 
average short lifetimes (~ 1 — 10 min) of the macrospicules. 
For example, the case shown in Fig. ([1} is characterized by 
the typical timescale of the self-heating, t ~ tA/Cs being of 
the order of 5mm. On the other hand, during one complete 
cycle of the self heating, almost 50% of the flow kinetic en- 
ergy will be extracted. For Fig. the same energy fraction 
is extracted from the mean flow, but the heating timescale 
would be ~ 2mm. Rather different situation is shown in 
Fig. (131, where the viscosity is not enough to saturate the 
instability and as is seen from the plot, the perturbed en- 
ergy is increased asymptotically. Nevertheless, the presence 
of viscous damping is still signiflcant, because as we see from 
Fig. ([3|, the self-heating rate for r = 20 equals ~ 500. On 
the other hand, we have already mentioned that ?7 ~ 1, thus 
the considered mode can extract the mean flow kinetic en- 
ergy in approximately 7mm. As is seen from these flgures, 
the more complex is the kinematics of the flow, the more 
efficient is the self-heating process. 

The results of the present study, being of a qualitative 
nature, indicate that nonmodal shear flow instabilities, cou- 
pled with the presence of the viscous dissipation, may en- 
sure the net chromospheric heating by the agency of the 
flow patterns of sufficiently complex kinematic structure. 
The presence of magnetic fields cannot alter this qualita- 
tive picture. On the contrary, it is known that in magne- 
tized helical flows nonmoda l shear instabilities t end to be- 
come much more powerful ijRoeava et al.ll2003al lbll: and it 
has been shown that se lf-heating in th e MHD may be quite 
efficient both generally (hi et al. 2006') and in the astrophys- 
ical context (Shcrgclaslivili ct al. 2006). It makes us believe 
that if self-heating is able to heat nonmagnetic part of the 
chromosphere of a solar-type star, it will certainly do the 
same job in a magnetized chromosphere too. We plan to 
study this issue in detail in the near future. 

In the standard, Biermann-Schwarz schild, scenario 
of heating o f the solar chromosphere (|Biermanl 1 19461 : 
ISchwarzschildl Il948l ) the velocity amplitude of upward- 
propagating acoustic waves grows due to the exponential 
decrease of the mass density with height; leading, in turn, 
to shock formation and dissipation of the acoustic waves 
and heating of the medium. The important feature of this 
mechanism is the exponential dependence of the mass den- 
sity on height, because for a medium with p — const this 
classic scenario does not work. Alternatively, the nonmodal 



self-heating works when the density is homogeneous! On the 
other hand, it crucially depends on the presence of the kine- 
matically complex shear flow. It is reasonable to expect that 
in realistic situations, viz. in chromospheric structures where 
both density stratiflcation and flows are present, these two 
mechanisms will be complementary and will serve as addi- 
tional sources of the heating. 

One of the goals of our future research will be also the 
study of the propagation of acoustic gravity waves in grav- 
itationally stratifled chromospheric flows in the presence of 
viscous dissipation. Both linear and nonlinear regi mes of 
acous tic-gravity wave propaga tion in two- (Kalkofc n et al.l 
|1994| ) and three-dimensional l|Bodo et all [2000 ) hydrody- 
namics and in the presen ce of the temperature inhomo- 
geneity (|Bodo et al.l 1200 ll ) have been previously studied. 
It is interesting to study the dynamics of these waves in 
the presence of velocity shear and viscous damping and to 
check whether the nonmodal self-heating influences these 
processes too. 

In order to study the role of the shear-induced non- 
modal instability for the problem of the chromospheric heat- 
ing in slowly rotating, solar-type stars we considered the 
model of a simple, quite general velocity pattern and wrote 
down hydrodynamic equations governing the evolution of 
acoustic perturbations within this flow. We have solved these 
equations for several, representative examples of flow pat- 
terns and found out that the instability and the viscous 
damping might lead to quite efficient self-heating of the 
"parent" flows. Analyzing the kinetic energy budget of a 
typical chromospheric macro-spicule, we found out that the 
self-heating might extract a significant part of kinetic energy 
from the velocity pattern in time-scales of the same order as 
the lifetimes of the corresponding structures observed in the 
solar chromosphere. Therefore, it is quite logical to believe 
that self-heating of chromospheric flow patterns contributes 
to the overall heating of the chromosphere by converting a 
considerable part of their kinetic energy into heat and ulti- 
mately destroying these flow patterns. 
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